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Pulsed neutron spectroscopy is used to make absolute measurements of the dynamic magnetic 
susceptibility of insulating YBa2Cu306.i5- Acoustic and optical modes, derived from in- and out-of-phase 
oscillation of spins in adjacent CuC<2 planes, dominate the spectra and are observed up to 250 meV. The 
optical modes appear first at 74 ± 5 meV. Linear- spin- wave theory gives an excellent description of the data 
and yields intra- and inter-layer exchange constants of 7y = 125±5meV and J± = ll±2meV respectively 
and a spin- wave intensity renormalization Z y = 0.4 ± 0.1. 



PACS numbers: 74.72.Bk, 61.12.-q, 74.25.Ha, 75.30.Ds 
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The fundamental building blocks of the high-T c superconductors are (nearly) square lattices of Cu 
atoms linked to their nearest Cu neighbours via the p-orbitals of intervening oxygen atoms. The blocks occur 
in units of one, two or three layers, which, in turn, are separated from each other by spacers serving as charge 
reservoirs. In general, the superconducting transition temperature rises with the number of Cu02 layers in 
the units. This means that multilayer materials, most notably YBa2Cu306+ x (T C =93K), have received much 
more attention than single-layer materials such as La2- x Sr x Cu04 (T c =38 K). Even so, the magnetic 
fluctuations, which are popular as a potential source of high-T c superconductivity, have been characterised 
over a much wider range of frequencies (co) and momenta (Q) for La2- x Sr x CuC>4 [1-5] than for 
YBa2Cu306+ x [6-9] . Indeed, while La2_ x Sr x CuC>4 , for both insulating and superconducting compositions, 
has been studied throughout the Brillouin zone for h(0 up to 0.4 eV, magnetic excitations in YBa2Cu3C>6+ x 
have only been observed for h(0 comparable to, or below, the pairing energy (40 meV) and thus much less 
than the likely coupling 7y responsible for the antiferromagnetism of the insulating parent. We have 

consequently performed magnetic neutron scattering experiments to establish the spin dynamics of 

YBa2Cu306+ x on a high energy scale. The sample chosen for this first measurement is an insulating 

antiferromagnet with ^ = 0.15. Our study reveals spin-wave excitations up to 0.25 eV. Below 
h(Og =74±5meV, we observe only acoustic spin waves where pairs of neighbouring spins in adjacent 

planes rotate in the same sense about their average direction. Above #G) g , we have discovered optical spin 

waves [10], where the spins in adjacent planes rotate in opposite directions, thus sensing the restoring force 
from the inter-planar coupling J± . Our data are well described by a linear spin-wave model in which an 
overall quantum renormalization of the spin- wave intensity Z^=0.4 ± 0.1 ( = 1 for the classical limit) is 

included and the intra-planar and inter-planar exchange constants are 7y = 125 ± 5meV and 7j_ = 1 1 ± 2meV 

respectively. 
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Experiments were performed on the HET spectrometer at the ISIS pulsed spallation source of the 
Rutherford-Appleton Laboratory. The experimental techniques for measuring high-frequency spin 
fluctuations using a direct-geometry chopper spectrometer such as HET are described elsewhere [3]. The 
sample used in the present investigation was a single crystal of YBa2Cu306.i5 with mass 96 g and room 

o o 

temperature lattice constants ofa = b = 3.857 A and c= 11.84 A corresponding to a Neel temperature of 400 

K [6]. In common with other large single crystals of YBa2Cu3C>6+x, our sample contains Y2BaCuC>5 as an 

impurity phase which constitutes up to 15% of the volume. We index reciprocal space using the tetragonal 
unit cell so that antiferromagnetic Bragg peaks occur at ^,0) and related positions, momentum transfers 

Qx'QyQz) m un its of A"l are then at reciprocal space positions (h,k,l) = [Q x aj2n, Q y b jln, Q z c/ 2 7i). 
Measurements were made with the (1 1 0) plane coincident with the horizontal scattering plane of the 
spectrometer. The detectors are 300-mm-long 25-mm-diameter tubes at 2.5 or 4 m from the sample with 
scattering angles in the range 2.5 < 20 < 28 deg. Absolute unit conversions were performed using a 
vanadium standard [10]. Data were collected at ambient temperature (T = 296 K). 

For the purpose of this experiment, we model YBa2Cu3C>6+x as a set of weakly coupled Cu02 
biplanes, i.e. a square-lattice bilayer antiferromagnet. In this case, the spin dynamics can be described using 
the Heisenberg Hamiltonian for a single bilayer 

# = X> s / .s J -+£y 1 s / .s / , (i) 

y y" 

where the first term represents the nearest-neighbour intra-planar coupling between Cu-II spins in the same 
plane and the second the nearest- neighbour inter-planar coupling between Cu spins in different planes. The 
introduction of the second term in (1) leads to two branches in the spin-wave dispersion [6] which are 
labelled according to whether nearest-neighbour spins in different planes rotate together (acoustic mode) or 
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in opposite directions (optical mode) about their average direction. Conventional linear-spin-wave theory of 



the Holstein-Primakoff type yields the imaginary part of the dynamic susceptibility (per formula unit) for the 



two modes as: 



i-y(G) 



i+r(Q)+j ± /2j n 
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2 2' 
COS 

J 



^AzQ z y(o)±0) op (Q)) (2) 



and 



n 



l-y(Q) + 7 1 /27 ll 

i+r(e) 



2 • 2 
sin 



(I 

v2' 



AzQ z }8{o)±(O ac (Q)), (3) 



where the dispersion relations are 

h(o(Q) = 2/,(l - Y 2 (2) + J±/2J\\ (1 ± y(Q))) 
(the plus sign is for the acoustic mode), y(Q) = l/2|cos(a<2 x ) + cos(a<2y)j and Az( =3.2 A) is the separation 

of the Cu02 bilayers. With the coupling terms in Eq. (1), there is no dispersion along the z-direction, only a 

modulation of the amplitude of the dynamic susceptibility, which may be used to distinguish the two modes. 
Further, the energy of the optical branch shows a minimum at Q = of M) g = 2- s JJ\\Jj_ . 

By varying the incident neutron energy Ej, we can arrange to cut across the spin waves at given 

energy transfers fi(0 and values of Q z (I). Figs. l(b)-(e) show data collected as a function of Q parallel to 

2 (1 \ 

(110) direction for various energy transfers h(0 such that sin —AzQ z ~1 i.e. / = Q z c/2n =1.8,5.5,.... 

V2 J 

Under these conditions, we are sensitive only to modes with acoustic character. At all frequencies, we 
observe a peak near Q={j^j) corresponding to the spin waves. Twin peaks, due to counter-propagating 

spin-wave branches, are not seen due to the poor out-of-plane resolution in the (l 1 0) direction. Figs. l(g)-(j) 



show similar data collected such that cos z 



(1 

v2' 



AzQ 7 «1 i.e. / = 3.7,7.3,. 



In this case, we are sensitive 



only to optical modes. Little or no magnetic scattering is evident below 82 meV. (The area under the 65 

meV peak is approximately 10 times smaller than for 82 meV). Our observations are in agreement with 
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of I 

thermal neutron scattering studies [6-8], made with 7zC0<42meV, which observed a sin —AzQ z 

modulation of the magnetic scattering at low energy transfers. In order to make a more quantitative 

interpretation we use our data to calculate the 2D local- or wavevector-integrated- susceptibility 
X2d(G z >G>) = JV'(2><fl) d 2 Qj^d 2 Q where the Q-integrals are over [Q x ,Q y ^ only and Q z is chosen to be 

close to an optical or acoustic position. The results are shown in Fig. 2. When the data are plotted in this 
way, it is clear that for energies #C0>82meV, the optical and acoustic intensities are equal within 

experimental error. On the other hand, for fi(0 < 65meV, the integrated response at the optical position is 
consistent with zero. Thus M) g lies between 65 and 82 meV. 

A more accurate estimate of the optical gap can be made from Fig. 3. This spectrum shows data 

collected as a function of ho) with Ej = 110 meV and integrated over wavevectors with 0.45 < h < 0.55, i.e. 
over the spin-wave cones emanating from Due to energy and momentum conservation in the 

scattering process, Q z varies with ha> in the way indicated by the horizontal axes of Fig. 3. Thus, Fig. 3 
displays the local susceptibilities associated with the acoustic and optical modes, weighted by their in- 
dependent structure factors. At lower frequencies (ha>< J\\), the 2-D local susceptibility associated with the 
acoustic mode is essentially co-independent and that associated with the optical mode has a step at h(O g (see 

solid lines in Fig. 2). Consequently, in Fig. 3 we expect to observe a simple sinusoidal modulation with Q z 

due to the acoustic mode, %2T>{Q Z ^) x sin 2 — Az <2 Z 1 for co<co g . For co>co g , a second sinusoidal 

v 2 J 

modulation of almost equal amplitude and ^ out of phase, due to the optical mode, is superposed, yielding a 

Q z - (and CO-) independent signal. By inspection, we can see that above approximately 70 meV the intensity 

is constant and below this value the characteristic acoustic modulation is observed. A resolution-corrected fit 
of Eqs. (2)-(4) to the data yields a value [10] for the optical gap of h(0„ = 74 ± 5 meV. 



The optical mode gap is simply proportional to the geometric mean of 7y and 7^ , and so by itself, its 
measurement establishes neither J± nor J\\. Indeed, to determine 7y and J^, it is necessary to determine not 
only ft(Og, but also the spin-wave dispersion as a function of in-plane momentum. We have therefore 

collected data at energy transfers sufficiently large to make the dispersion obvious. Fig. 4 shows the results 

obtained for an incident energy Ei = 600 meV. For fi(0=\15 ± 15 meV, the scattering, though broadened by 
the dispersion, is still peaked at {^,\ ,/)• For %(0=2A5 ±15 meV, the dispersion causes the scattering to be 

peaked near (0.35,0.35,/). Finally, for slightly larger energy transfer ha> = 275 ± 15 meV, there is a much 
diminished modulation, consistent with a cut-off of approximately 250 meV for single-magnon scattering. 
The solid lines in Figs. 4(b)-(e) show resolution-corrected fits of linear spin-wave theory. From fitting each 
energy individually, we determine the dispersion relations shown in Fig. 4(a). A simultaneous fit including 
the lower frequency data (see solid lines in Figs. (l)-(3)) yields exchange constants 7y = 125±5meV and 

J ± =ll±2meV. 

A range of values for J\\ have been inferred from other measurements. Rossat-Mignod and co- 
workers [7] deconvolved reactor-based measurements of low-frequency spin waves and found a spin-wave 
velocity, near of c=1000 ± 50 meVA, which if we assume the absence of further neighbour 

interactions (confirmed by the present experiment), implies 7y ~ 200 meV. Shamoto and co-workers [8] 

o 

performed a similar study and found c = 655 ± 100 meVA. This is in agreement with the initial slope of the 
acoustic dispersion curve we have measured. The only pre-existing information on the magnetic energy 
scale at short wavelengths is from interpretation of the two-magnon Raman effect [11], which suggests, for a 
model including only a single near neighbour exchange, 7y =120 meV. Our value for 7y is substantially 

below that (0.16 eV) [3] for La2Cu04. This discrepancy may be due to greater off- stoic hiometry in our 

YBa2Cu306.i5 sample as well as its larger lattice constants. 
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It is well known that in two-dimensional quantum antiferro magnets, the spin-wave energies and the 
overall intensity of the spin waves show quantum renormalizations Z c and respectively [12,13] when 

compared with their classical (large S) values. Neutron scattering cannot directly measure Z c , thus the bare 

exchange constants must be inferred from our effective or measured coupling constants using theoretical 
estimates [12] of Z c ~1.18. In contrast, we are able to measure Z^ experimentally by placing our 

measurements on an absolute intensity scale. Following this procedure, we obtain a value of Z^=0.4 ±0.1, 

close to that given by a | spin-wave expansion [13] and indistinguishable from that given by our 

measurements [5] on La2Cu04. 

Not as much is known about J± as about 7y and Z^. Even so, detailed analysis [14] of the 

essentially zero-frequency magnetic response measured [15] by nuclear resonance for Y2Ba4.Cu70i5 

suggests 5<7^<20 meV. It is also interesting to ask what happens to J± in the metallic and superconducting 

samples, especially because one explanation [16] for the celebrated spin gap behaviour of YBa2Cu306+ x , 

which is not so clearly seen in the single layer compounds, requires J^^O. Magnetic scattering for all 

superconducting samples of YBa2Cu306+ x displays the modulation characteristic of acoustic spin-wave 
modes [7,9]. Because scattering has been observed up to h(0 ~ 45 meV = 0.6/zC0 g (x=0.15), the bilayer 

coupling for even the optimally doped superconductor is unlikely to be much smaller than 4 meV or 45 K 
which is 36% (=0.6 2 ) of J ± (x = 0.15) . 

In summary, we have measured, for the first time, single-magnon excitations in the insulating 
antiferromagnet YBa2Cu3C>6.i5 throughout the Brillouin zone. Our results are well described by a linear 
spin-wave model for a 2D square bilayer, provided an appropriate overall amplitude renormalization is 
included. The resulting inter- and intra- planar exchange constants are 7|| = 125±5meV and 

7 = 11 ± 2meV respectively. 
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Figure Captions 

FIG. 1. (a) and (f) show the dispersion relations for acoustic and optical spin-wave branches respectively. 
Remaining panels are energy scans showing magnetic scattering from YBa2Cu3C>6.i5 for wavevector 
transfers along Q = (hhl). (b)-(e) / is chosen to give scattering from acoustic modes, (g)-(j) / is chosen to 

give scattering from optical modes. Incident energies and Q z wavevector components are Ei = 600, 1 10, 75, 
75, 600, 110, 75, 75 meV and 1=5.3, 5.6, 5.5, 1.9, 3.6, 7.2, 7.1, 3.7 for (b)-(e), (g)-(j) respectively. 

FIG. 2. The 2-D local susceptibility as a function of energy transfer obtained by integrating over spin-wave 
peaks and correcting for the Cu 2+ magnetic form factor, Bose factor and instrumental resolution, (a) acoustic 
positions and (b) optical positions. The figure is a compilation based on the analysis of data such as those in 
Fig. 1 with various incident energies. 

FIG. 3. Data collected with Ej = 110 meV, integrated over the spin-wave cones in the (<2 x >Qy) plane near 

Q = (^-,^,/)and corrected for the Cu 2+ magnetic form factor. The resulting spectrum is collected along a 

trajectory in (Q z ,h(o) defined on the upper and lower axes. Under these conditions the sinusoidal intensity 
modulation with Q z will disappear at %(£> g ~lA meV. Solid line corresponds to the best description in terms 

of spin wave theory. Dotted line is the acoustic mode contribution. 

FIG. 4. (a) The dispersion relation obtained from independent fits at each energy transfer (all data in paper). 
Closed circles and solid line are for acoustic modes, open circles and dotted line are optical modes, (b)-(e) 
Constant energy scans showing high-frequency magnetic scattering from YBa2Cu3C>6.i5. Data were 
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collected with ki II (001) and Ej = 600 meV. Counting time was 29 h at 170 (J,A proton current with a Ta 
target. Solid lines are resolution-corrected fits of a linear spin-wave model for a bilayer (see text). 
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